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Abstract 

We study motion in the field of two fixed centres described by a family of 
Einstein-dilaton-Maxwell theories. Transitions between regular and chaotic 
motion are observed as the dilaton coupling is varied. 
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The classical two centre problem describes the motion of a small mass in the field of two 
fixed centres. The solution for motion restricted to a plane was given by Euler in 1760, and 
the general solution was found by Jacobi in 1842. In the Newtonian case the centres 
may be kept fixed by balancing gravitational attraction against electrostatic repulsion. The 
relativistic analog of the two centre geometry was found independently by Majumdar 
and Papapetrou and was later shown Q to describe two or more extremal Reissner- 
Nordstrom black holes. The general relativistic two (or A^) centre spacetime has now been 
extended to form a one-parameter family of solutions 0]. These solutions contain a dilaton 
field (j) in addition to the gravitational field g and Maxwell field F = dA and are described 
by the action 

S = j (fx [R - 2{d(t)f - e-^^-^F^) . (1.1) 
The static A^-centre solution solution takes the form 

A = , dt , 

e-^ = H^/^^^"'^ , (1.2) 
where if is a harmonic function describing the positions Xj of the masses Mf. 

(1 + a'^)Mi , , 

i ' 

Each black hole has mass Mj, electric charge Qi and dilaton charge Sj. These charges satisfy 
the extremality condition = M'^ + T? . The parameter a labels the family of solutions, and 
is related to various reductions of iV = 8 supergravity. Some interesting special cases are: 
a = Einstein-Maxwell; a = 1/a/3 Einstein-Maxwell reduced from 5 — >^ 4 dimensions; a = 1 
string theory; a = Einstein gravity reduced from 5 — 4 dimensions . Solutions with 
a = 0, 1/v^, 1, arise in type II string theory as marginally bound states of elementary 
solutions with a = v^- Much of the recent interest in these solutions has focused on the 
duality between extremal black holes and intersecting D-branes [Q. This broader context 
is not the focus of our current study. Rather, we are interested in using ( |1.2| ) to study the 
relativistic two centre problem. 

The Newtonian two centre problem leads to equations of motion that are separable 
and hence integrable |Q. Chandrasekhar considered the Einstein-Maxwell two centre 
problem but was unable to integrate the equations of motion. Methods borrowed from 
dynamical systems theory were then used to prove that null and timelike geodesies of the 
Majumdar-Papapetrou spacetime were chaotic Here we extend these results to include 
null and timelike geodesies of the general Einstein-Maxwell-dilaton two centre spacetimes. 
Most of the paper is devoted to null geodesies as the chaotic dynamics admit a complete, 
analytic description. The null geodesies can be thought of as describing ultra-relativistic 
chaotic scattering. Timelike geodesies and the motion of extremally charged test particles 
are discussed at the end of the paper. Extremal test particles have charges that satisfy 
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= im? + 0"^, where m is the mass, q is the electric charge and a is the dilaton charge. 
The motion of an extremally charged test particle in the field of two massive fixed centres 
provides a first approximation to motion in moduli space - the space of all static solutions. 

It is well known that null geodesies in static spacetimes correspond to ordinary geodesies 
of a three dimensional optical metric which, in our case, is given by 

dsl = H^/^'-^'^'Ux' (1.4) 

The optical metric is defined on M^/{xj} and it is complete if < 1. If < 1 we have 
an outer asymptotically fiat region connected to a number of asymptotically fiat (a = 0) or 
asymptotically conical (0 < a < 1) regions surrounding the centres. The asymptotic regions 
are separated by Einstein-Rosen type 'throats'. If = 1 these throats become infinitely 
long tubes. If > 1 then the optical metric is incomplete and the centres are singularities 
of the optical metric at finite optical distance. 

Rather remarkably, we find that null geodesies in the Kaluza-Klein limit a = ^/S are 
non-chaotic, and can be integrated using methods familiar to Jacobi For intermediate 
values of a we find that the chaotic bands in phase space grow as a increases from a = 0, 
reaching a peak at the string value a = 1, and then shrinking to zero as a ^ y/S. 

The story for uncharged (g = cr = 0) timelike geodesies is not so rich as chaos reigns 
for all values of a. In contrast, the character of the dynamics for extremal test particles 
varies strongly with a. Extremal test particles behave very much like photons. A moving 
extremal particle never comes to rest. When moving, extremal particles interact via purely 
velocity dependent forces - there is no static component to the attraction. The transition 
pattern between regular and chaotic motion mirrors that of null geodesies. In particular, we 
find that the motion of a small extremal black hole in the field of two larger black holes is 
integrable in the Kaluza-Klein limit. 

II. REGULAR AND CHAOTIC MOTION 

A. Time and chaos 

It was the best of times, it was the worst of times. 
Dickens, The tale of two cities 

When studying the dynamics of generally covariant theories, one no longer has the rigid 
Newtonian concepts of absolute space and universal time. Time and space become relative 
concepts, and standard measures of chaos that rely on the metrical properties of phase space 
become observer dependent. In our "tale of two centres" there is no choice of time coordinate 
that is any better or any worse than any other. 

For this reason, we abandon gauge dependent measures such as Lyapunov exponents and 
metric entropy in favour of invariant measures such as fractal dimensions and topological 
entropies. 
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B. Fractal methods 



The geodesies will be chaotic if there exists a chaotic invariant set of orbits. This set is 
usually referred to as a strange repellor or strange saddle. The terminology "strange" refers 
to the set having some non-integer multi-fractal dimensions. The term "repellor" or "saddle" 
indicates that the orbits are unstable in some directions and stable in others. Uncovering 
these fractal structures in phase space provides a gauge invariant way of showing that the 
motion is chaotic [p|,p!0| 

Of particular interest is the future invariant set. For unbound motion, the future invariant 
set correspond to those trajectories that approach the two centres from infinity with an 
impact parameter that allows them to take up everlasting periodic orbits. Since these 
orbits never exit the scattering region, they cannot be assigned to a particular asymptotic 
outcome. Thus, the future invariant set forms the boundary between the various outcomes. 
In a numerical experiment the future invariant set can be uncovered by studying the nature 
of the boundaries between the outcome basins. 

The experimental setup is shown schematically in Fig. 1. The two centres, both with 
mass M = 1, are placed at x = (0, ±1, 0). Since angular momentum is conserved about the 
y-axis, it is enough to consider motion restricted to the plane z = 0. Null geodesies are fired 
into the scattering region from the point (xci/o) at an angle 6q from the x-axis. The initial 
velocities are given by 

u^ = uy = (2 1) 

The trajectories are evolved numerically until an outcome is reached. We assign four 
different outcomes on the basis of where the trajectory ends up in the (x, y) plane: 

1(0,1), 11(0,-1), III(±oo,oo), IV(±oo,-oo). (2.2) 

The first two of these outcomes will not occur for all values of a, since the capture cross 
section of each centre goes to zero for a > 1. Initial conditions were chosen by setting xq = 10 
and selecting [yo, Oq) from a 500^ grid. Points in this grid are colour coded according to their 
final outcome using the colour scheme: I dark grey; II light grey; III black; IV white. The 
results of the numerical experiment are displayed in Fig. 2 for three values of a. The first 



two graphs show evidence of "chaotic gravitational lensing" |Tl|. 

The boundaries between the outcomes appear to be fractal for a = 0, 1 and smooth for 
a = ^/S. To confirm these impressions, portions of each image are enlarged in Fig. 3. Again, 
the boundaries for a = 0, 1 are clearly fractal, while the boundaries for a = -\/3 look quite 
smooth. 

As a further check, the capacity or box counting dimension of the images is computed 
using the method described in Ref. We find 

D{a = 0) = 1.36 ±0.02, 
D{a = 1) = 1.50 ±0.02, 

L)(a = v^S) = 1.01 ± 0.02 . (2.3) 
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Within numerical tolerances, these dimensions support the contention that null geodesies in 
geometries with a = 0, 1 are chaotic, while the a = ^/S geometry has regular null geodesies. 
Repeating this analysis for several values of a we find that geometries with a < y/S dis- 
play chaotic scattering, while those with a > \/3 do not. The chaotic behaviour is most 
pronounced at the string value a = 1. 



C. Curvature methods 



As an alternative to our fractal methods, several groups p2|-p!5|] have advocated curvature 
as a coordinate independent tool for forecasting chaos. The idea is to extend Hadamard's 
classic result that the geodesic fiow on a compact manifold with all sectional curvatures 
negative at every point is chaotic. So far, these attempts have failed to yield any reliable 
method to forecast chaos. The various criteria proposed are neither necessary nor sufficient 
for predicting chaos. 

Despite these shortcomings, the idea of using curvature methods is not entirely without 
merit. With a little care it is possible to arrive at a necessary, but not sufficient, criteria for 
chaos. The reasoning is as follows: For chaos to occur the phase space must contain a chaotic 
invariant set. Since the elements of this set are unstable periodic orbits, the dynamics must 
admit such trajectories. If there are no unstable periodic orbits, then the dynamics will 
not be chaotic. Applying this test to geodesic motion requires us to show (1) the manifold 
admits periodic geodesies, and (2), most of these orbits are unstable. A sufficient criteria for 
condition (2) to hold can be given in terms of orbit-averaged sectional curvatures. Chaotic 
behaviour is ruled out if either (1) or (2) is not satisfied. 

To improve our test so that it is both a necessary and sufficient for chaos would require 
some notion of mixing. For Hadamard, the mixing comes about because the hyperbolic 
manifold is compact. In general, the geodesic flow will not be restricted to a compact 
region, so we have to look for other mixing mechanisms. For near-integrable systems the 
mixing can be caused by a homoclinic or hetroclinic tangle, the existence of which can be 
probed using the Melnikov method 

Condition (1), i.e. the existence of periodic null geodesies, is equivalent in a static space- 
time to the existence of closed geodesies in the optical metric. Some partial information 
is provided by the Benci-Giannoni theorem The Benci-Giannoni theorem guarantees 



the existence of at least one non-constant closed geodesic of a complete (but not necessarily 
compact) riemannian manifold provided a condition on the fall-off of the sectional curvatures 
and a topological condition hold. The optical metric is complete if a < 1 and the topological 
condition holds in that case. The sectional curvature condition for the optical metric is easily 
seen to hold if a < 1. Thus, perhaps surprizingly, the Benci-Giannoni theorem guarantees 
the existence of at least one, and probably very many, exactly periodic null geodesies no 
matter how many centres one has and no matter how one positions them. 

Condition (2) can be checked using the geodesic deviation equation. Returning to the 
usual spacetime metric, the geodesic deviation equation describes how nearby geodesies 
separate: 

= -R^^Xn'W^ . (2.4) 
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We demand that the deviation is a spacehke four vector orthogonal to the four velocity 
u^, ie. 



n^u^ = M% = 0, -1 n% > 0. (2.5) 

It is convenient to introduce a set of non- rotating orthonormal basis vectors, {a = 0..3), 
with Eq set equal to u" [Q. The remaining three basis vectors can then be used to describe 
the spacelike deviation vector n'^ = rfE^ where i = 1..3. Only two basis vectors are required 
to describe the deviation when the four velocity is null [^]. The advantage of this approach 
is that it separates the rotation of the geodesic congruence from its spreading. This allows us 
to write the deviation equation ( |2.4| ) in terms of ordinary, rather than covariant derivatives: 

^ = -R'o^on^ ■ (2-6) 
Contracting (p.6| ) with rij and averaging over an orbit we find 

^ ' ' = -< K{u, n > +^ 2.7 

where n = (n'^n^Y^'^, 

<K{u,n)>= ^— - — , (2.8) 



and 



</>= <Pfds= l' fds, A/ = /° 



(2.9) 

The quantity < K{u, n) > is essentially the sectional curvature in the plane spanned by u'^ 
and n^, averaged over one orbit. If <K{u, n) > is negative we have 

A{dn'^/ds) > 0. (2.10) 

That is, the rate of deviation increases with each orbit: 

A periodic orbit is unstable if any of its three (two for null geodesies) average sectional 
curvatures is negative. 

To see that conditions (1) and (2) only provide a necessary condition, we can apply the 
test to a single black hole spacetime with a = 0. In this case there is a circular photon orbit 
at r = 2M (in areal rather than isotropic coordinates) with four velocity 

= (1,0,0,—). (2.11) 

4A/ 

A suitable pseudo-orthonormal set of basis vectors are 

Eo^ = n^ E3^ = (4,0,-,0) 

M 

E^ = (0, i, 0, 0), E^ = (0, 0, — , 0) . (2.12) 

2M 
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Both Eq and E3 arc null and satisfy q^^EqE^ — —1. Since wc arc dealing with null geodesies, 
we need only consider deviations in the plane spanned by Ei and E2. Writing the deviation 
vectors as 1 = /Ei and q = we find 

The rate of contraction in the 1 direction exceeds the rate of expansion in the q direction 
due to Ricci focusing. This focusing is due to the term 

in the Raychaudhuri equation for the expansion, 9 — Despite the focusing term, the 
circular photon orbit is unstable against radial perturbations. Thus, a single extremal 
Reissner-Nordstrom black hole spacetime satisfies the necessary conditions for chaos to occur. 
However, geodesies are integrable in this spacetime, so our curvature condition is not a 
sufficient condition for chaos. 

Returning to the geodesic equation for Z'*, we have 

Here we have chosen the affine parameter for the null geodesic to coincide with the coordinate 
time t. Solving for l{t) we find 

l^kexpf^^] . (2.16) 



One might be tempted to say that the radial direction has a positive Lyapunov exponent A = 
l/(-\/32M), but this statement is highly coordinate dependent. For example, an observer 
free falling into the black hole sees the trajectories diverge at a rate /(r) ~ r^^, so she 
would conclude that A = 0. To avoid this type of ambiguity, we define unstable to mean 
there is at least one negative orbit-averaged sectional curvature. 



III. SYMBOLIC DYNAMICS 



Taken on their own, the numerical results provide a solid, but rather unenlightening 
description of the dynamics. A far more satisfying description can be found using symbolic 
dynamics. Symbolic dynamics describes the topology of trajectories in phase space. Because 
the detailed local dynamics does not enter into this description, the symbolic dynamics can 
be studied analytically, even when the trajectories are not integrable. 

Since we are considering chaotic scattering, the unstable periodic orbits that form the 
strange repellor are of particular importance. Typical scattering trajectories experience 
chaotic transients as they pass through the scattering region. The dynamics of these tran- 
sients is completely encoded by the strange repellor. Thus, by uncovering the symbolic 
coding of the repellor, we learn a great deal about the dynamics. 
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To find the coding we place windows in phase space, positioned in such a way that 
topologically distinct periodic orbits pass through the windows in a distinct order. Each 
time an orbit passes through a window, the symbol for that window is recorded. The 
resulting string of symbols is the symbolic coding of the orbit. A unique symbolic coding 
can be found by demanding that each distinct physical orbit has a distinct symbolic coding, 
and that every symbolic coding describes a physical orbit. 

Once the coding has been found, the symbolic complexity of the repellor can be measured. 
If the dynamics is chaotic, then the number of periodic orbits will grow exponentially as the 
symbolic length of the orbits is increased. 



A. Allowed Orbits 

The primary closed null geodesies of the two centre spacetime are shown in Fig. 4. Orbits 
of type a) and b) require that each centre is capable of bending a trajectory through at least 
Ab = 2ti. Similarly, type c) requires A;, > vr, and type d) requires Af, > vr. 

Since large angle scattering occurs for trajectories with small impact parameter 6, ie. 
close to one of the two centres, the maximum scattering angle produced by each centre can 
be approximated by ignoring the distant centre and considering a one-centre spacetime. The 
spherical symmetry of the one-centre spacetime reduces the dynamics to one dimensional 
motion in an effective potential V{r): 

where h is the impact parameter and 

V{r)=(l+^-^-^^] r-\ (3.2) 



r 

The scattering angles can be calculated from the equation 



The result is 



/■Umax 2h dll 

Al= / ^ -vr. (3.4) 

Jo V(l + (1 + a2)MM)«/(i+'^') - uW 

Here u = 1/r and "Umax denotes the point of closest approach. The integral ( |3.4|) will be 
finite unless the denominator admits a double root. A double root occurs when 

2M(1 + (1 + a2)Mu)(i-'^')/(i+"') = h , (3.5) 

and 

(l + (l + a2)MM)2/(i+"') =m6. (3.6) 
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Solving for u and h we find 



^--=m[t^2)^ (3-7) 

and 

6erit = 2M(^— -^j . (3.8) 

It is no coincidence that these are the same values of u and h for which give rise to unstable 
photon orbits. We see that unstable photon orbits are only possible if a < 1. Thus, the 
scattering angle is finite for all a > 1. The limiting case a = 1 has to be handled 
separately. A direct evaluation yields 



-4M2 



1 



2b , f2M\ , , 

sin^i M- • (3-9) 



- 4M2 



Here the critical impact parameter separating capture from scattering is 6crit = 2M. If we 
write h = 2M(1 + e^) where e <^ 1, then we find 

~ 7r-4. (3.10) 

e 

Hence the scattering angle can be infinite when a = 1. This means a glory is possible even 
though the effective potential does not have a turning point. 

For all a > 1 the scattering angle will be finite. For a just a little larger than 1 the 
scattering angle can still be large enough to allow several temporary photon orbits. This 
partial glory adds to the complexity of the allowed periodic orbits, but not to the same 
degree as a full glory. 

Unstable photon orbits will still be possible in the two centre spacetimes so long as 
Afo > 71. Note that the largest scattering angles occurs for orbits that pass very close to 
one of the centres. For these orbits our approximation is especially good. We find that the 
critical value of A;, is reached when a = ^/S. For small b the scattering angle reaches 

Af = 7c- — + ^--.... (3.11) 
* M 12M3 ^ ^ 

In this case there are can be just one unstable periodic orbit in the two centre spacetime, and 
thus no chaos. Later we will show that the Kaluza-Klein limit a = ^/S is actually integrable 
for null geodesies. 

The existence of periodic null geodesies and the associated glories in these black hole 
spacetimes should be constrasted with their complete absence in cosmic string spacetimes 
|T9| . The reason for this difference is that for cosmic strings the sectional curvatures are 



non- negative |T9|. This fact again points to the importance of the sectional curvatures in 
general relativity as a possible diagnostic tool for chaos or its absence. 
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B. The symbolic coding 



Using the primary unstable orbits as a guide, we see that it is natural to place three 
windows along the axis connecting the two centres. The placement and labelling of the 
windows is shown on the left of Fig. 5. Using these windows, all orbits can be represented 
by strings of O's and ±l's, with the restriction that no symbol follows itself. The allowed 
transitions are shown schematically on the right of Fig. 5. 

Since a complete orbit must contain an even number of symbols, the counting of or- 
bits is made easier if we shift to the new symbolic alphabet A = {(1, — 1), (— 1, 1)}, 
B = {(0, — 1), (— 1, 0)}, C = {(0, 1), (1, 0)}. We define the length of a symbohc sequence 
to be the sum of the number of incontractible loops around each centre. This number is a 
topological invariant as motion is restricted to the (x, y) plane. The labelling of the windows 
was chosen so that the length of an orbit is given by the sum of the absolute values of the 
symbols used to describe the orbit. For example, the primary orbit d) has the symbolic 
coding 1, 0, —1, = CB, and is thus of length k = 2. Here overline denotes a sequence to be 
repeated. 

The counting of orbits as a function of their length is a simple exercise in combinatorics. 
An orbit of length k = 2p + n will be made up of n + p symbols {A, B,C}, where p is the 
number of A's, i the number of i?'s and n — i the number of C"s. Thus, the number of orbits 
at order k is given by 



[k/2] n=k~2p . . 



p=0 i=0 

+ + (3.12) 

Here the notation \x\ denotes the integer part and the double sum is over a trinomial 
combinatoric factor. Readers familiar with number theory will recognise + 1 as the silver 
mean 2 + 1/(2 + 1/(2 + .... The expression for Nilz) looks more natural when derived from 
a recurrence relation. Longer orbits can be generated from shorter orbits by inserting A's, 
-B's and C's, so that 

N{k^2)=N{k)^2N{k^\) . (3.13) 

The first term comes from inserting an A, and the second from inserting a i? or a C. A 
direct counting of orbits reveals A^(l) = 2 and iV(2) = 5. Using these initial values, one 
immediately arrives at ( 3.12| ) as the solution to the recurrence relation ( p.l3 ). 



C. Topological entropy 

Since the symbolic coding is based on an uneven three symbol coding (the A's are twice 
as long as the i?'s and C's), we expect the topological entropy to lie between that of a 
straight two symbol coding and a straight three symbol coding, ze., In2 < Ht < In 3. This 
is indeed the case: 
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Ht = lim y In N{k) =\n{V2 + l) . (3.14) 

Unlike the metric or Kolmogorov- Sinai entropy, the topological entropy provides a coordinate 
independent measure of chaos in general relativity |[T0|j2H| . 

For a > 1 the symbolic coding starts to get pruned as there can only be a finite number 
of orbits around each centre between excursion to the second centre. This limits the number 
of 5's or C's that can be strung together. If the maximum scattering angle about one centre 
is Amax, then the longest string of the form 5" or C" is given by n = [Amax/7i"]- Once A^ax 
drops below 2n, no orbits of type a) or b) can be inserted between orbits of type c) or d). 
The symbolic coding can then be reduced to a binary alphabet and the topological entropy 
drops to Ht = In 2. The topological entropy then remains unchanged until a = a/3, at 
which point A^ax ^ and orbits of type d) are no longer possible. The symbolic coding 
collapses to a single letter and the topological entropy drops to zero. When the topological 
entropy of the strange repellor vanishes, the dynamical system is non-chaotic. This analysis 
is in complete agreement with the numerical results of the previous section. 

D. The integrable limit a = \/3 

We have shown that the scattering of massless particles in the Kaluza-Klein two centre 
geometry is not chaotic. The repellor has integer capacity dimension and zero topological 
entropy. Now we will show that the Kaluza-Klein two centre problem is integrable for null 
geodesies. 

Null geodesies are best studied in the optical metric 

da^ = -de + H^/^^+'''\dx^ + dy^ + dz^) . (3.15) 

Following Jacobi |l[] we adopt prolate spheroidal coordinates 

X = sinh ip sin 6 cos 0, y = sinh ip sin 6 sin 0, 

z = cosh ip cos 9 . (3.16) 

In these coordinates, the harmonic function H is given by 

H =1 + W/Q, Q = sinh^ tfj + sin^ 6, 

W ={Mi + Ms) cosh + (Ml - Ms) cos^ . (3.17) 



The Hamilton- Jacobi equation. 



takes the form 



{dt) Hm+a^)smh^^|JsmH{^<j)) if4/(i+a2)Q U 9^ J + ' ' 1 0.(3.19) 
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Since t and are cyclic coordinates, their canonically conjugate momenta, E and L, are con- 
stants of the motion. The Hamilton- Jacobi equation will be separable, and hence integrable 
if 

S = ~Et + + S^ij) + Se{9) . (3.20) 

Substituting this anzatz into ( |3.19| ) we find the system only separates if a = y/S. The 
integrable limit a = a/3 is characterised by an addition constant of the motion a: 

dS \^ r2 

^ =« -^ + S'(sinh2V' + (Mi + M2)cosh7/'), (3.21) 

dip ) smh ■?/' 

and 



a^y ^ il 

09 " sin^ 9 



-a - + E\sm^9 + (Mi - Ms) cosh^) . (3.22) 



IV. TIMELIKE TRAJECTORIES 

Test particle with mass m, electric charge q and dilaton charge a obey the equation of 
motion 

^ + Viyu^ = ^e-"<^F>^ - —(u^u^Vy + V°)</) . (4.1) 

aA ^'^ m ^ m 

A curious situation occurs for extremal test particles initially at rest. They are characterised 
by 



q = Vl + a^m, a = am, u' = H^'^^+''\ (4.2) 

and 



= _/72/(i+a ^ VlT^HEi - aV'<t) . (4.3) 
CiA 



Inserting the background solution (|1.2|) into (^]3|) yield; 



s 

Initially static extremal test particles remain at rest. On reflection this is not so surprising. 
An extremal test particle can be thought of as a small extremal black hole moving in the 
field of other, larger black holes. When at rest, the small black hole acts like another centre 
in the static multi-black hole solution. 
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A. Uncharged test particles 



The simplest timelike trajectories are the timelike geodesies followed by uncharged (g = 
(T = 0) test particles. Depending on their energy, uncharged test particles may be confined 
to the region of space near the two centres. Those that are not confined will eventually be 
captured by a black hole or scattered to infinity. Using the same techniques we applied to 
null geodesies, it is easy to show that unconfined timelike geodesies are chaotic for all values 
of a in the range < a < \/3. 

For a single centre the motion is described by 

(^^y^E^.Virf, (4.5) 
where the effective potential has the form 

V = THtL^ (h'"'*''' + ^) ■ (4.6) 



Here E and L are the test particle's conserved energy and angular momentum per unit 
rest mass. Trajectories with E > 1 are able to escape to infinity. As we found for the null 
geodesies, massive particles with non-zero angular momentum can only be captured if a < 1. 
Moreover, there are no unstable periodic orbits if a > 1, but there are stable orbits for all a. 
For two centres, separability of the Hamilton- Jacobi equation, 

dS dS o , , „s 

q — TT-j = , (4.7) 

is broken in the Kaluza-Klein case (a = v^) by the term 



H^l'' . (4.8) 



As the energy of the test particle is increased, the non-integrable term ( |4.8| ) can be neglected 
so that ultra-relativistic test particles have non-chaotic geodesies. 

An interesting new ingredient enters into the motion of confined test particles. There will 
be a locus of points surrounding the two centres where confined geodesies will momentarily 
come to rest before falling back toward the centres. This locus of points is called the zero 
velocity curve. In a two centre spacetime the zero velocity curve makes possible a whole 
new range of unstable periodic orbits not possible for scattering trajectories. The new class 
of unstable orbits comprises all possible traverses between points on the zero velocity curve. 
As a result, confined geodesies are very much more chaotic than unconfined geodesies. 

The highly chaotic nature of confined geodesies is dramatically illustrated in Fig. 6. 
The basin structure results from treating timelike geodesies in a Kaluza-Klein two centre 



spacetime as a Hamiltonian exit system Since a > 1 and < 1, the geodesies can 

neither escape to infinity nor fall into a black hole. In other words, there are no asymptotic 
outcomes. To remedy this problem, we introduce exits around both centres. A small circle 
was drawn around each centre. When a trajectory exits the system an outcome is assigned 
according to which centre the exit enclosed - black for (0, 1), white for (0, —1). Test particles 
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were released from rest with initial positions in the (x, y) plane taken from a 500^ grid. The 
initial grid point was colour coded according to its outcome. 

It might appear that we are forcing a square peg into a round hole by using fractal basin 
boundaries in a situation where there are no natural asymptotic outcomes. It is certainly 
true that a system with bound trajectories is well suited to study by standard techniques 
such as Poincare surfaces of section. A Poincare section would also provide coordinate 
independent, fractal information about the dynamics, but without the need for exits to be 
introduced. What we wanted to show was that fractal basin boundaries work for both bound 
and unbound systems. On the other hand, Poincare sections can only be used for bound 
orbits and this means they are of very limited use in general relativity |10| . 



B. Extremal test particles 

Earlier we showed that static extremal test particles remain at rest. The converse is also 
true: If an extremal test particle is in motion, it will never come to rest. (We are neglecting 
retardation effects due to the emission of gravitational, scalar or electromagnetic radiation). 
In loose terms, extremal test particles only feel velocity dependent forces. 

In many ways, an extremal test particle in motion behaves like a photon. In fact, it can 
be shown that extremal test particles moving in the metric ( |1.2| ) follow null geodesies 



of the five dimensional metric 



The five dimensional Hamilton- Jacobi equation can be reduced to the four dimensional form 



2| 

The dynamics of extremal test particles in the various two centre geometries can be 
charted using the same techniques we applied to null geodesies. A numerical survey in- 
dicated that the pattern of dynamical behaviour exhibited by extremal test particles is 
essentially identical to what we found for null geodesies. This correspondence appeared to 
be independent of velocity, so long as the particles were moving. For example, the outcomes 
basins of Fig. 7 should be compared to those displayed in Fig. 2. The morphology of the 
basins is identical. 

Once again we see that the Kaluza-Klein geometry appears to admit regular trajectories. 
Adopting prolate spheroidal coordinates, it is a simple exercise to show that the Hamilton- 
Jacobi equation ( [4.10| ) is separable when a = \/3- The only extra terms in the separation 
equations (|3.21|) and ( p.22| ) are 2mE sinh^ and 2mE sin^ 6 respectively. 
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V. CONCLUSIONS 



Using a combination of fractal and topological techniques, we have given a complete 
description of the Einstein-dilaton-Maxwell family of two centre problems. Unlike their 
Newtonian counterpart, most relativistic two centre problems are not integrable. There are 
arguments indicating that for more than two centres the Newtonian problem is not Liouville 
integrable [24| and it therefore seems likely that in our case as well the motion will not be 
integrable for more than two centres. 

We did find one exceptional case where our tale of two centres came full circle to its 
Euler-Jacobi antecedent: the motion of massless particles and small extremal black holes 
is integrable in the field of two fixed centres residing in a Kaluza-Klein compactified five 
dimensional spacetime. Thus, we can add the Kaluza-Klein two centre problem to our 
meagre collection of integrable three body problems. 
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FIGURES 




FIG. 1. Initial conditions for the null geodesies. 
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FIG. 5. The symbolic dynamics. 
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FIG. 6. Outcomes for massive, uncharged test particles released from rest in the (a = two 
centre geometry. 
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